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Space-Time Curvature Signatures in Bose-Einstein Condensates 
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We derive a generalized Gross-Pitaevski (GP) equation immersed on a electromagnetic and a 
weak gravitational field starting from the covariant Complex Klein-Gordon field in a curved space- 
time. We compare it with the GP equation where the gravitational field is added by hand as an 
external potential. We show that there is a small difference of order gz/c? between them that could 
be measured in the future using Bose-Einstein Condensates (BEC). This represents the next order 
correction to the Newtonian gravity in a curved space-time. 
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Among the four forces of nature, the gravitational 
force is the hardest one to study. It is orders of mag¬ 
nitude weaker than the other three forces and requires 
astronomical masses to see corrections beyond the ba¬ 
sic Newtonian formula. The lab scale masses that we 
can manipulate introduce gravitational forces that are 
very hard to measure accurately. This is reflected in 
the poor precision that we currently have on the de¬ 
termination of the Newtonian constant of gravitation 
(G = 6.67384(80) x lO^kg^nrV 2 ) Q. 

Most of the experiments that measure the gravita¬ 
tional force use basically the same method pioneered by 
Cavendish 0. Atomic interferometers have been recently 
introduced as an alternative way to implement sensitive 
gravimeters [3]. Here, cold atoms are launched vertically 
and subjected to a sequence of laser pulses that split, re¬ 
combine and interfere the atoms. Using a Bose-Einstein 
Condensate (BEC) reduces the ballistic expansion of the 
atomic cloud and allows for a larger fall time thus in¬ 
creasing the precision of the gravimeters. 

Sensitive atomic gravimeters are useful in the detec¬ 
tion of gravitational waves fH, tests to Einstein’s general 
relativity , study of short range gravitational forces [|j 
and coherent evolution of delocalized quantum objects 
[?]. Measurements of g (the gravitational acceleration of 
the earth) typically reach 9 digits of precision. This has 
been recently improved by two order of magnitude by 
launching a BEC in free fall in a 10 m tower [8j. Planned 
improvements to this include using a 146 m tower |9j, 
sending the atoms in parabolic flights M or putting the 
atoms in space EH- 

BECs are usually described by a Gross-Pitaevskii (GP) 
equation, where an external field is added in order to con¬ 
fine the system or to include the earth’s gravitational ac¬ 
celeration ([U|)- Nevertheless, the GP equation is non- 
relativistic and non-covariant; it is not invariant under 
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Lorenz transformations unlike the electromagnetic field. 
A GP equation obtained this way could only show New¬ 
tonian effects of gravity. A more complete description is 
obtained by treating the gravitational field not as a force, 
but as the curvature of space-time as in Einstein’s theory 
of relativity. In Ref. [13l | it was shown that the Klein- 
Gordon equation in a curved space takes the GP form 
after a simple transformation and includes extra terms 
that are interpreted as the gravitational field and a fi¬ 
nite temperature contribution. The gravitational com¬ 
ponents include terms beyond the Newtonian gravity. In 
this work we propose that these terms could become the 
focus of future experimental measurements. 


I. KLEIN-GORDON IN A WEAK 
GRAVITATIONAL FIELD 

References [l4| and [l5| derive the relativistic theory 
of a BEC in a quantum field theory description. Starting 
from the Klein-Gordon equation in a flat space-time, a 
generalized GP equation is obtained for relativistic and 
finite temperature fields. The generalization for an ex¬ 
panding universe is given in [13]. Here we apply the same 
approach on a curved space for a BEC living in a weak 
gravitational field, the Klein Gordon equation for a com¬ 
plex scalar field is given by 


(i) 

where $ is the complex scalar field and A^ is the cor¬ 
responding electromagnetic four vector. In this work we 
use the Mexican hat scalar field potential given by 

m 2 r 2 A 

U=^$* + _|$| 4 . (2) 

We add an external interaction U ex t that can represent 
for example the potential that confines the condensate 
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like a laser, a box or the gravitational field seen as an 
external one. The metric of a weak gravitational field is 

M 


ds 2 = — (1 + 2il>)dt 2 + (1 — 2cj))gijdx l dxi , (3) 

where ip, (p are the gravitational potentials and gtj is 
the 3-dimensional flat-space metric. The Klein-Gordon 
equation with this metric reads 

n NEW $ - +ie<&\7nA^ - e 2 A 2 $ + 2ieA fJ Vu$ = 

d<P* 

a NEW® - ^-j—■ + *e[(l + 20)V ■ A — -(1 — 20)0 

+ V(0 — (p) ■ A + -(ip + 30)y>]$ 
c 

- e 2 ((l + 20)A 2 - (1 - 2iP)ip 2 )<i> 

+ 2ie((l + 2(p)A ■ V<F - -(1 - 20)^4>) 

c 

Orri ^ 

= -J^U^ (4) 

with e = e/(Hc ), A and <p the magnetic vector and elec¬ 
tric scalar potentials such that A M = (ip, A). The metric 
modifies the D’Alambertian as 

g 2 

a NEW = -(1 — 20) c2 Q t2 + (1 + 20)V 2 

+ ^[0 + 3 </,]A+V(^-</»)-V (5) 

The metric not only produces the gravitational potential, 
but it also affects the kinetic energy term, something that 
does not happen in the traditional treatment that adds 
the gravitational potential by hand. 

In a curved space-time in the weak field limit the 
Maxwell equations read 


V ■ E + V (0 — 30) ■ E = p 
E + (0 - 30)E + B x V(0 - 30) - V x B = J 

V ■ B - 4B ■ V0 = 0 
B - (1 + 20)V x E - 2E x V0 = 0 (6) 

where p and J are the charge density and electric current 
respectively and the electric (E) and magnetic (B) fields 
are defined by the relations 


(see also G2) 


h 2 A 

i(i — 2 <p)hi> = ——^new^ + 


me 2 {(p + Uext)^ - 


imr(( 1 + 20)V ■ A-(1 — 21 p)ip + 

c 


2m 


V(0 — (p) ■ A + - {ip + 3 <p)ip) — 
c 


2m 


((1 + 2(p)A 2 — (1 — 2ip)ip 2 ) 


T - 


i—[( 1 + 2(P)A • V4- - (1 - 2ip)<p(-V - ™$)]. (8) 

TO C 


where to = mc/h. This last equation is again the Klein- 
Gordon equation (Eq. [XI), but written in terms of the 
function T. 

If we remove in Eq. [8] the gravitational fields <p and 
ip and we set the charge e = 0, we recover the Gross 
Pitaevskii equation in the non-relativistic limit. There¬ 
fore we interpret Eq. [8] as the generalization of the GP 
equation for a relativistic charged bose particle in elec¬ 
tromagnetic media in a weak gravitational field. 

In Einstein’s theory with a weak gravitational field we 
have <p = ip. In some alternative theories of gravity this 
identity does not follow fl8i |. For simplicity, we assume 
a static gravitational field (<9 t 0 = 0 = 0) like that of the 
Earth. The result after taking the Newtonian limit is 


i(l-20)M' = 

fc2 \ 

-(1 + 20)-V 2, P H-I'FI 2 4' + mc 2 04' + 

2 m 2m 

mc 2 U ext ’F - i^r(( 1 + 20)V • A - -(1 - 20)0)4' + 
2m c 

^~r((l + 20)A 2 — (1 — 20)</? 2 )4' - 
i^[( 1 + 20)A • VIE' - (1 - 2 <p)p{^A> - irh&)\. (9) 

We interpret Eq. [9] as the generalization of the GP 
equation in a gravitational field immersed in a electro¬ 
magnetic potential. 

We first study the difference between the traditional 
approach to Newtonian gravity to the curved space-time 
derivation used in this work. For that we neglect the 
electromagnetic field (e = 0). Taking 0 = 0, we recover 
the flat space result that corresponds to the usual GP 
equation 


(1 - 20)2 

1 (d A \ 

~ (l-20)(l + 20) V 7 

To write down the Klein-Gordon equation (Eq. 0j) in its 
GP form we apply the transformation $ = 4'e _ * mc t ! n 


h 2 A 

iM> = - — V 2 + + —14'| 2 4' + mc 2 U ext '& (10) 
2m 2m 

Taking U ex t = gz/c 2 we obtain the gravitational po¬ 
tential of the Earth ( mgz ). Here the gravitational in¬ 
teraction is added as an external potential. This case 
represents the Newtoninan version of a BEC in a gravi¬ 
tational field. 
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In the version of Einstein we introduce the gravita¬ 
tional potential through the space-time curvature. Tak¬ 
ing U ex t = 0 and 0 = gz/c 2 gives the gravitational poten¬ 
tial as well, but there is a correction since 0 appears also 
in the term of Eq. [9] To determine the size of the cor¬ 
rection lets consider an atomic interferometer. Here an 
atom (or the BEC wavefunction) is splitted in two, and 
the two components evolve in free fall for some time be¬ 
fore recombining them. The separation is usually smaller 
than 1 cm |§], but it is conceivable to achieve separations 
up to 1 m in the near future by combining large free fall 
times with large momentum transfer. A z =1 m separa¬ 
tion gives a 0 ~ 10 -16 , that is small compared to the 1 
in the d' term of Eq. [9] 

The combination gz/c 2 is the perturbative parameter, 
in analogy to the term v 2 / c 2 of special relativity. The 
observation of the correction requires doing atomic inter¬ 
ferometry with 16 digits of precision. The current record 
is at 11 digits [I], but the projected sensitivity of the 
space project is at 15 digits of precision 0- Having the 
atoms in space goes in the wrong direction since that 
eliminates g. Instead, the present work suggests mov¬ 
ing towards atomic interferometry in strong gravitational 
fields. 

It is possible to reach higher values of g by taking ad¬ 
vantage of the equivalence principle. Instead of the grav¬ 
itational field of the Earth one could accelerate the com¬ 
plete experimental setup. This is clearly complicated, 
but in the case of acceleration introduced by Bloch oscil¬ 
lations things might be simpler. Here it is only necessary 
to sweep the phase of the two counter propagating lasers 
beams to produce an accelerated lattice [191 ]. There is a 
price to pay in the precision of the measurement due the 
reduced measurement time since, for example, a linear 
acceleration of 100 g gives already a displacement of 5 nr 
in 0.1 s. An alternative is to have the BEC in a centrifugal 
force. Here the gravitational field is (/ = v 2 /c 2 ln(R/R 0 ), 
with v c the tangential velocity of the BEC, R the radius 
of the centrifuge and Rq an arbitrary reference radius. 
For R = 1 m and v c > 3.2 m/s the centrifugal force is 
bigger than the gravitational force of the Earth. 

We apply a Madelung transformation 'k = yJnc lS to 
Eq[9]with e = 0, with n representing the density number 
of particles and S the velocity super potential 

v=—VS. (11) 

m 

The real and imaginary parts of Eq. [9] in terms of the n 
and S variables are 


(1 - 2 <f>)h + (1 + 20)V • (nv) - (1 - 2 <j>)j = 0(12) 

r 2 


(1 — 2 < f >)- + -(1 + 20 ) — + ( 

c 2 c A 


+ - 


A 


h 2 □ 


NEW 


y/n 


2m 2 2 m 2 c 2 yfn 

respectively, with the flux j = nv and 


1 V 

2 (1 ^c 


= 0 ( 13 ) 


J = n- 


tS = n- 


(14) 


We interpret Eq. [T2]as the generalized continuity equa¬ 
tion in gravitational and electromagnetic fields. This 
equation differs from the one derived using pure fluid 
mechanics in the factors in front of the density and ve¬ 
locity terms. Equations fl2l and fl3l are the Klein-Gordon 
equation in a weak gravitational field, in other words, 
they are the Einstein-Klein-Gordon equations written in 
the variables n and S. We calculate the gradient of Eq. 
m to obtain the corresponding momentum equation 

(1 — 20)nv + (1 + 20)ra(v ■ V)v = 

2 

nF ext + n(l + 2— \-2—^-)F ( / > — 
c c 

—Vp + nF q + Vcr (15) 

where p is the pressure with a state equation p = 
(A /m 2 )n 2 , Fq = — V((/i 2 /277i 2 )V 2 v / n/v(n) is the quan¬ 
tum force and a is the viscosity (see [2l| for details of 
the non-gravitational case). As expected, the external 
force (by unit of mass) is defined by F ext = —S7U ex t and 
the gravitational force by = — c 2 V0. Equations [T21 
and fl5l correspond to the Klein-Gordon equation (Eq. [D 
written in terms of n and v. 

The difference between the two procedures to include 
the gravitational potential becomes evident from Eq. [15] 
If we consider the gravitational field as an external force 
(non-covariant case), then space-time is flat, 0 = 0 and 
thus F^ = 0. In this case the hydrodynamic equa¬ 
tion corresponds to the traditional one. If instead we 
use the covariant form of the equations we must set 
F ext = 0 and include the gravitational force in the F^ 
term. The coefficient in front of both terms differs by 
2v/c+ 2 v 2 / c 2 ~ 2 v/c. The modification is very small 
because of the v dependence. To connect with the pre¬ 
vious treatment, the wavefunction phase S evolves at a 
different rate at different heights z because of the grav¬ 
itational potential, so that S = mgzt/h. Using Eq. [Uj] 
we have 2 v/c = 2 gz/c 2 which is the same parameter we 
obtained before. 

In summary, we present a field theoretical approach to 
describe a BEC in a curved space. The derivation results 
in a generalized GP equation in a gravitational and elec- 
tromagneitic fields. We identify the expansion parameter 
gz/c 2 that gives corrections that could be verified with 
precise enough atomic interferometry. This makes BECs 
in strong gravitational (or accelerated) fields a very inter¬ 
esting candidate to study corrections beyond Newtonian 
gravity. 
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